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ABSTRACT 


The  complex  permittivity  ( s )  and  permeability  ( ju  )  of  a  material  detennine  the 
response  of  the  material  to  electromagnetic  radiation.  Usually,  the  real  parts  of  s  and  // 
are  positive  for  naturally  occurring  materials  at  microwave  frequencies.  Metamaterials 
are  engineered  media  that  are  designed  to  have  either  a  negative  permittivity  or 
penneability  or  both.  Negative  permeability  and  negative  permittivity  would  cause 
electromagnetic  waves  traveling  through  this  medium  to  exhibit  unusual  characteristics 
such  as  power  flow  in  a  direction  opposite  to  the  phase  velocity.  In  this  thesis,  the  wave 
matrix  approach  is  used  to  calculate  the  total  reflection  and  transmission  coefficients  of  a 
multilayered  structure.  The  method  is  applicable  to  all  types  of  materials,  including 
metamaterials.  Several  layered  configurations  are  studied  including  both  metamaterial 
and  conventional  dielectric  layers.  A  MATLAB  program  is  developed  to  examine  the 
effects  of  frequency,  angle  of  incidence  and  polarization.  The  results  are  compared  to 
published  data.  Potential  applications  of  metamaterials  are  also  discussed. 
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I.  INTRODUCTION 


A.  BACKGROUND 

For  frequency  domain  analysis,  all  materials  can  be  described  electrically  by  their 
complex  permittivity  ( s  )  and  permeability  ( // ).  These  parameters  detennine  the  response 
of  the  material  to  electromagnetic  (EM)  radiation.  Usually,  the  real  parts  of  s  and  //  are 
both  positive  in  naturally  occurring  materials.  While  some  materials  could  have  negative 
real  s  (such  as  conductors  at  optical  frequencies),  no  natural  material  which  has  a 
negative  real  p  is  known.  Over  30  years  ago,  the  Russian  physicist  Victor  Veselago 
predicted  the  existence  of  electromagnetic  materials  called  metamaterials,  which  have 
negative  pennittivity  and  permeability. 

Metamaterials  are  engineered  media  that  are  designed  to  have  either  a  negative 
permittivity  or  permeability  (single  negative  material,  SNG)  or  both  (double  negative 
material,  DNG).  They  are  created  by  using  composite  structures  that  are  generally 
periodic  with  periods  much  less  than  a  wavelength.  Negative  penneability  and  negative 
pennittivity  would  cause  electromagnetic  waves  traveling  through  this  medium  to  exhibit 
certain  unusual  characteristics.  For  example,  waves  are  refracted  at  a  negative  angle 
relative  to  conventional  materials.  Also,  due  to  the  negative  index  of  refraction,  the  phase 
and  group  velocity  vectors  of  an  electromagnetic  wave  are  in  opposite  directions.  These 
unusual  properties  may  be  used  to  build  new  devices,  like  highly  sensitive  antennas  or 
even  a  perfect  planar  lens. 

In  this  study,  the  propagation  characteristics  of  planar  multilayered  structures 
(including  both  metamaterial  and  conventional  dielectric  layers)  are  examined  using  a 
wave  matrix  approach.  The  reflection  and  transmission  coefficients  are  calculated,  and 
applications  are  discussed  for  anti-reflection  and  high-reflection  structures.  The  effects  of 
frequency,  angle  of  incidence  and  polarization  are  examined. 

B.  OBJECTIVES 

The  objective  of  this  thesis  is  to  develop  a  matrix  approach  to  predict  the 
scattering  characteristics  of  planar  multilayered  structures  composed  of  metamaterial, 
dielectric  and  magnetic  layers.  The  wave  matrix  formulation  is  developed  to  calculate  the 
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reflection  and  transmission  coefficients  of  multiple  layers  of  arbitrary  thicknesses  and 
materials.  Therefore,  with  the  help  of  MATLAB  software,  the  effects  of  frequency,  angle 
of  incidence,  polarization  and  layer  configuration  can  be  examined.  To  validate  the 
MATLAB  program,  many  test  cases  are  presented  and  simulation  results  are  compared  to 
published  data. 

C.  ORGANIZATION  OF  THESIS 

In  Chapter  II,  equations  for  transmission  and  reflection  coefficients  for  a  single 
boundary  between  two  dielectric  media  are  presented.  The  wave  matrix  approach  for 
layered  media  is  described.  The  reflection  and  transmission  coefficient  equations  are 
given  for  layered  media. 

In  Chapter  III,  metamaterials  and  their  unusual  characteristics  are  discussed.  The 
potential  applications  of  DNG  materials  are  presented. 

In  Chapter  IV,  the  propagation  characteristics  of  multilayered  structures 
consisting  of  dielectric  slabs,  metamaterials,  or  both,  are  discussed.  Different  methods  are 
described  to  achieve  zero  reflection  or  high  reflection.  The  effects  of  frequency,  angle  of 
incidence  and  polarization  are  shown.  The  potential  applications  of  these  structures  are 
described. 

Finally  in  Chapter  V,  the  results,  conclusions  and  suggestions  for  future  studies 
are  discussed. 
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II.  WAVE  REFLECTION  AND  TRANSMISSION 


A.  INTRODUCTION 

Electric  charges  generate  electric  fields  and  electric  currents  generate  magnetic 
fields.  If  the  charge  and  currents  sources  were  to  vary  with  time,  the  electric  and 
magnetic  fields  become  unified,  and  the  link  between  them  generates  electromagnetic 
waves  (EM)  that  can  travel  through  free  space  and  in  material  media.  If  a  plane  wave 
from  one  medium  meets  a  different  medium,  it  is  partly  reflected  and  partly  transmitted 
depending  on  the  constitutive  parameters  of  the  two  media  involved.  This  chapter  begins 
with  an  examination  of  the  reflection  and  transmission  properties  of  plane  waves  incident 
on  planar  boundaries,  and  continues  with  the  development  of  wave  matrices  for  layered 
media. 

B.  REFLECTION  AND  TRANSMISSION  COEFFICIENTS 

Consider  a  unifonn  plane  wave  traveling  along  the  +z-direction  that  is  normally 
incident  at  the  boundary  between  medium  1  (/ux,£x)  and  medium  2  ( ju2,£2 ).  The  incident 
(i),  reflected  (r),  and  transmitted  (/)  waves  are  shown  in  Figure  1. 


x 

A 


Medium  #1  (Ju]  ,£,  ) 


Medium  #2  ( jil2 ,  s2  ) 


(Incident  wave) 


a 


(Transmitted  wave) 


(Reflected  wave) 


(z  =  0) 


Figure  1 .  A  plane  wave  incident  normally  on  an  interface  between  two  different 

media  (After  Ref.  [1]). 
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The  incident  wave  is  traveling  in  the  +z  direction  in  medium  1.  For  nonnal 
incidence,  the  incident  electric  field  vector  can  be  expressed  as 


E  ,0)  =  E  ,0er*ax 

and  the  magnetic  field  vector 

H,(z)  =  H,„e^a>,=5»c-»-.=S 

7, 

where  //,  is  intrinsic  impedance  of  material  1 ,  given  as 


(2.1) 


(2.2) 


(2.3) 


and  p  =  //'  -  jff ,  e  =sr'  -  js"  are  the  relative  permeability  and  permittivity  of  the 
medium.  Also,  yx  is  the  propagation  constant  of  the  medium,  and  is  given  as 

y1=al+  .  (2.4) 

The  reflected  wave  is  traveling  in  the  -z  direction  in  medium  1.  The  reflected 
electric  field  vector  can  be  written  as 

Er(z)  =  Er0e^dx  (2.5) 

and  the  magnetic  field  vector  is 

H r(z)  =  Hr0ey,z(-d  )  =  -%He**ay .  (2.6) 

7, 

The  transmitted  wave  is  traveling  in  the  +z  direction  in  medium  2,  and  its  electric  field 
vector  is 

Et(z)  =  Et0e-^dx  (2.7) 


and  its  magnetic  field  vector 
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(2.8) 


H  ,(z)  =  H  l0e~r2=av  =  —e^av 
Vo 

where  rj2  and  y2  are  defined  in  manner  similar  to  Equations  (2.3)  and  (2.4). 

In  Figure  1,  the  total  field  in  medium  1  includes  both  the  incident  and  reflected 
fields.  In  medium  2,  only  the  transmitted  field  is  present.  Then,  the  total  fields  can  be 
written  as 


E>=E,+Er 

Hi=H,+H„ 

E2=E? 

h2=h(. 


The  planar  boundary  located  at  z=  0  separates  the  two  media.  At  the  boundary,  the 
tangential  components  of  the  electric  and  magnetic  fields  are  continuous.  Thus, 


Hi,  -  H2,  -»  H;0  +  Hr0  -  H,0  -» 


i  -  Ef0 

(2.9) 

(E?o  —  Er0)  _  E/0 

(2.10) 

7I\ 


V  2 


From  Equations  (2.9)  and  (2.10), 


E 


r  0 


(2.11) 


E 


to 


2*72 
V2  +77i 


E,0. 


(2.12) 


Now  it  is  possible  to  obtain  the  reflection 
coefficient  ( r  )  from  Equations  (2.11)  and  (2.12), 

r_E,0  _7l2~7l\ 
E,o  rl2  +rl\ 

.  _  E,o  2tj2 

E/0  r?2+rll 


coefficient  (T)  and  transmission 


(2.13) 

(2.14) 
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To  simplify  the  analysis,  normal  incidence  is  assumed.  Situations  that  are  more  general 
are  considered  in  the  following  chapters. 

C.  WAVE  MATRICES  FOR  LAYERED  MEDIA 

Thus  far,  the  focus  has  been  on  uniform  plane  waves  traveling  through  two  media 
and  a  single  boundary,  next,  a  more  complex  situation  is  considered;  that  is,  multilayered 
media.  First,  the  two-media  equations  are  rewritten  in  matrix  form  and  then  the  matrix  is 
extended  for  N  layers. 

Figure  2  shows  incident  waves  approaching  the  boundary  from  both  sides 
(traveling  in  the  +z  and  -z  directions),  and  waves  traveling  away  from  the  boundary  in 
both  media.  By  denoting  the  +z  traveling  incident  wave  bx ,  and  -z  traveling  incident 


incident  wave  ( b] ) 

Transmitted  wave 

~\ 

Reflected  wave 

r  ◄ - 

r  ^2 

Reflected  wave 

r  J 

1  |  Transmitted  wave 

- ►  J 

Incident  wave  ( r, ) 

h 

Figure  2.  Normal  incidence  at  a  boundary  between  two  media, 
wave  r2 ,  then  incident  waves  become, 

ri=TA+T2ir2  (2.  i  5) 

and 

b2=r2r2+rl2bl  (2.16) 
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where  the  first  subscript  on  r  denotes  the  incident  medium  and  the  second  subscript 
denotes  the  transmitted  medium.  Now  it  is  possible  to  rearrange  Equations  (2.15)  and 
(2.16),  and  write  them  in  the  matrix  form  [Ref.  2], 


bl 

1 

f  1 

_r  ^ 

1  2 

b2 

_rl_ 

hi 

vri 

712721  _  ^1^2  y 

_  ri 

(2.17) 


The  wave  matrix  in  Equation  (2.17)  is  for  a  two-layered  structure.  Next,  it  will  be 
extended  to  N  layers.  Figure  3  shows  forward  (+z)  and  backward  (-z)  traveling  waves 
through  the  N  layers. 
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Figure  3.  Forward  and  backward  traveling  waves  through  N  layers  (After  Ref.  [2]). 


To  simplify  the  equations,  the  forward  traveling  waves  are  named  bx,b2,...,bN  and 
the  backward  traveling  waves  are  named  rrr2,...,rN .  At  z  =  tx,  the  forward  traveling 
wave  becomes  bxe~Jnt' ,  and  backward  traveling  wave  becomes  rleJr,t' ,  where  yn  is  the 
propagation  constant  for  layer  n  defined  in  Equation  (2.4) 

co  I -  2n  I -  ,  . 

Yn  =  —y]Mr^rn  (2‘  18) 

where  X  is  the  wavelength  in  the  free  space. 
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For  the  plane  wave  in  Figure  3,  traveling  from  free  space  through  the  /Y-laycred 
medium,  the  wave  matrices  can  be  cascaded  to  give 


N  1 

=nf 


f  y„t„ 
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n=l  l'ii 


Te 


y„<„ 


-y„‘„ 


JN+ 1 
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\a2l 


11  M12 
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22  J 


JN+ 1 


(2.19) 


If  the  thickness  of  the  last  layer  of  the  medium  extends  to  oo ,  then  it  is  possible  to  state 
that  bN+l  =  0 .  The  total  transmission  and  reflection  coefficients  of  the  medium  become, 


■w+i 


_rL_a2L 

T  i 

bi  au 


(2.20) 

(2.21) 


Thus  far,  the  equations  in  this  chapter  are  based  on  nonnal  incidence.  The  wave 
matrix  can  be  extended  to  non-normal  angles  by  using  the  wave  impedance  in  place  of 
the  intrinsic  impedance.  The  nonnalized  wave  impedance  for  a  parallel-polarized  wave 
incident  at  an  angle  6*  from  the  interface  normal  is  given  as  [3] 


z_Vw-sin2fl 

sr  cos  6i 

and  for  a  perpendicularly-polarized  wave 

v_  fur  cosffi 

4Ert*r  “Sin2  6, 


(2.22) 


(2.23) 


The  transmission  and  reflection  coefficients  at  the  interface  between  the  two  media  are 
given  as 


7  -7 

Y-  2 

(2.24) 

Z 2  +  Zj 

r  =  l  +  r 

(2.25) 

where  Zx  is  the  incidence  medium  and  Z2  is  the  transmitted  medium.  Finally,  the  phase 
of  the  layer  (thickness  t )  is  given  as 
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(2.26) 


D.  SUMMARY 

The  first  part  of  this  chapter  presented  the  equations  for  transmission  and 
reflection  coefficients.  These  equations  are  based  on  nonnal  incidence,  and  a  single 
boundary  between  two  dielectric  media. 

The  second  part  described  wave  matrices  for  layered  media.  Then,  reflection  and 
transmission  coefficient  equations  are  given,  both  for  parallel-polarized  and 
perpendicularly-polarized  waves  incident  at  an  angle  6i .  These  situations  are  simulated 

with  MATLAB  6.5  in  Chapter  IV,  both  for  dielectric  media  and  metamaterials.  However, 
first  a  brief  description  of  metamaterials  is  given. 
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III.  METAMATERIALS 


A.  INTRODUCTION 

In  1968,  Victor  Veselago  (a  Russian  physicist)  forecasted  the  existence  of  a 
material  in  which  both  permittivity  and  permeability  were  assumed  to  have  negative  real 
values  (negative  index  of  refraction)  [4].  These  materials  have  been  called  metamaterials, 
backward- wave  materials  and  left-handed  (LH)  materials.  Although  these  materials  were 
not  known  to  exist  at  that  time,  he  concluded  that  with  these  negative-index  materials  all 
known  wave  propagation  and  optical  behaviors  would  be  changed. 

Materials  with  negative  permittivity  occur  naturally  (e.g.,  metals  at  optical 
frequencies),  but  those  with  negative  permeability  do  not.  Victor  Veselago’s  prediction 
was  verified  after  more  than  30  years  when  in  1999,  John  Pendry  (Imperial  College, 
London)  showed  how  these  materials  could  be  created  artificially  [4],  Then,  in  the 
following  year,  he  published  a  paper  in  which  he  proposed  that  metamaterials  could  be 
used  to  make  a  perfect  lens  [5]. 

B.  NEGATIVE  PERMITTIVITY  AND  PERMEABILITY 

Metamaterials  are  created  by  using  composite  structures  which  are  engineered  to 
have  properties  that  do  not  occur  naturally.  Figure  4  shows  one  type  of  negative-index 
material  that  was  constructed  by  David  Smith  at  the  University  of  California  in  2000  [4]. 


Figure  4.  An  artificial  negative-index  material  made  from  grids  of  rings  and  wires 

(From  Ref.  [4]). 
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The  cells  are  small  compared  to  the  wavelength  so  that  macroscopically  the  effective 
index  of  refraction  is  negative. 

Metamaterials  are  engineered  media  that  are  designed  to  have  either  a  negative 
pennittivity  or  permeability  (SNG)  or  both  (DNG).  For  a  DNG  material,  the  index  of 
refraction  for  real  permittivity  and  permeability  is 

n  =  -yl\/ir\\er\ ,  ^  <°’G  <°  (3.1 ) 

and  the  intrinsic  impedance 

7= Wkl/kl  •  (3-2) 

They  have  several  unusual  properties  when  compared  to  the  regular  materials  [6]: 

1 .  The  phase  velocity  of  EM  waves  is  anti-parallel  to  the  group  velocity  and 
the  flow  of  the  energy.  Thus,  propagation  is  described  by  a  left-hand  rule. 

2.  Waves  are  refracted  at  a  negative  angle  relative  to  conventional  materials. 

3.  Evanescent  waves  grow  with  distance  into  the  medium. 

4.  The  Doppler  shift  is  reversed  and  Cerenkov  radiation  propagates 
backwardly. 

Artificial  materials  are  a  mix  of  regular  materials  combined  to  obtain  specific 
characteristics.  A  negative  s'r  can  be  obtained  by  a  three  dimensional  array  of  wires  [7] 
(Figure  5).  The  effective  dielectric  constant  is  given  as 


Figure  5.  Three-dimensional  grid  of  thin  wires  approximates  a  plasma  (From  Ref. 

[7]). 
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where  co  =  — - - - — ,  r  is  the  radius  of  the  wire,  a  is  the  grid  spacing  and  cr  is  the  wire 

a  In (a/r) 

conductivity. 

A  coplanar  ring  (CPR)  is  an  example  of  a  configuration  that  influences  both  the 
effective  pennittivity  and  penneability  of  a  material  (Figure  6).  The  effective 
penneability  of  the  CPR  is  given  as 

_  2  /  2 

,  7tr  a 


corju0  7r2co2ju0Cr 3 


£  ( Is  | 

where  p  is  the  resistivity  of  the  metal  (ohms/m),  C  =  —  In  —  is  the  capacitance  per 

n  \t  ) 

meter  of  two  parallel  strips,  a  =  lattice  spacing  in  the  plane  of  rings  and  1=  spacing 
between  sheets  of  rings. 


Metal  rings 
with  gaps 


+  +  +T 


+  +  + 


Figure  6.  Coplanar  ring  (From  Ref.  [7]). 
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1.  Left  Hand  Rule 

All  materials  can  be  described  electrically  by  an  equivalent  complex  permittivity 
and  permeability.  The  permittivity  of  a  material  ( s  =  s0sr )  determines  its  response  to  an 

electric  field,  while  penneability  (//  =  //„//, )  determines  its  response  to  a  magnetic  field. 

Together  these  two  parameters  determine  the  reaction  of  the  material  to  the  EM  radiation. 
All  natural  (passive)  materials  follow  a  right-hand  (RH)  rule  because  of  their  positive 
values  of  relative  pennittivity  and  permeability.  However,  the  double  negative-index 

materials  follow  the  left-hand  (LH)  rule  [6].  Figure  7  demonstrates  these  rules  where  S  is 
the  poynting  vector  defined  as 

S=ExH* .  (3.5) 

A  material  for  which  the  RH  rule  holds  is  called  right-handed  material  (RHM).  The  LH 
rule  holds  in  a  left-handed  material  (LHM).  In  a  LHM,  the  phase  of  the  waves  (given  by 
vector  ak  )  moves  in  the  direction  opposite  to  the  direction  of  the  energy  flow  (given  by 

vector  S ) . 


s  >  0,//  >  0 


E 

Right-Handed  Materials 


Left-Handed  Materials 


Figure  7.  Illustration  of  vector  relationships  for  right-hand  and  left-hand  rules. 

2.  Negative  Refraction 

Snell’s  law  gives  the  relationship  between  angles  of  incidence  and  refraction  at  an 
interface  between  two  media  with  different  indices  of  refraction 

sin  6>,  =  ft2  sin  02 . 
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(3.6) 


The  same  law  applies  for  negative-index  materials,  but  because  of  the  negative  sign,  the 
refracted  wave  remains  on  the  same  side  of  the  normal  as  the  incident  wave.  In  Figure  8, 
a  wave  that  strikes  an  interface  between  air  and  a  positive-index  material  is  refracted 
towards  the  surface  normal,  while  in  a  negative  index  material,  it  is  refracted  on  the 
opposite  side. 


Air 

n>0 

oy 

Refracted  wave 

Incident  wave 

n<0 

0-V 

Refracted  wave 

Incident  wave 

Figure  8.  Refraction  in  a  positive-index  material  (top),  and  in  a  negative-index 

material  (bottom). 


Measurements  and  simulations  with  artificial  negative-index  materials  have  been 
performed.  In  Figure  9,  an  experiment  at  MTI  has  shown  the  predicted  refraction 
behavior  [8].  The  index  of  refraction  of  Teflon  and  a  LHM  are  determined  at  10.5  GHz 
by  measuring  the  angle  of  refraction.  The  propagation  is  toward  the  top  of  the  figure  and 
the  beam  is  normally  incident  on  the  prisms.  The  opposite  faces  of  prisms  are  cut  at  an 
angle  of  18° .  To  prevent  leakage  out  of  the  sides,  absorbing  foam  is  used  in  the  lower  half 
of  the  figures.  The  measured  indices  were  n  =  1.52  ±0.07  for  Teflon  (actual  index  for 
Teflon,  n  =  1 .5  )  with  a  refraction  angle  of  6  =  28°  ±  2°  and  n  =  -0.36  +  0.13.  For  the  left- 
handed  material,  the  measured  refraction  angle  was  0  =  -6.4°  ±  2.4° .  A  microwave  beam 
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traveling  upwards  is  refracted  in  the  positive  direction  when  it  leaves  a  Teflon  wedge, 
while  it  is  refracted  in  the  opposite  direction  when  it  leaves  a  metamaterial  prism 
constructed  from  wires  and  rings.  Two  lobes  appear  in  the  LHM  prism  due  to  the 
corrugation  of  the  refractive  surface. 


?o 


t 

incident  wa' 


8  12 
position  (cm) 


1C 


Figure  9.  A  microwave  beam  refracted  in  Teflon  (a),  and  refracted  in  a  metamaterial 

(b)  (From  Ref.  [8]). 


In  theory,  a  negative  refraction  material  can  be  used  to  build  a  completely  flat  lens 
that  focuses  light  to  a  point  (Pendry’s  perfect  lens  [5]).  Instead  of  making  specific 
concave  or  convex  surfaces,  metamaterials  should  be  able  to  serve  as  a  perfect  lens  with  a 
simple  plane  surface.  Unlike  standard  lenses,  metamaterial  lenses  have  no  optical  axis. 
Figure  10  demonstrates  the  idea  of  a  planar  metamaterial  lens. 
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Air  Metamaterial  Air 


Figure  10.  Planar  metamaterial  lens,  known  as  Pendry’s  perfect  lens  (After  Ref.  [5]). 


3.  Growing  Evanescent  Waves 

Evanescent  waves  are  formed  when  sinusoidal  waves  are  internally  reflected  off 
an  interface  at  an  angle  greater  than  the  critical  angle.  In  a  RHM,  the  evanescent  waves 
decay  exponentially  in  amplitude  with  distance  from  the  interface  at  which  they  are 
formed.  Left-handed  materials  can  cancel  the  decay  of  evanescent  waves  [5].  Pendry 
showed  that  evanescent  waves  emerge  from  the  far  side  of  the  medium  enhanced  in 
amplitude  by  the  transmission  process  [5].  This  does  not  violate  energy  conservation 
because  evanescent  waves  do  not  carry  energy. 

4.  Reversed  Doppler  and  Cerenkov  Radiation 

The  Doppler  Effect  is  the  apparent  change  in  frequency  or  wavelength  of  a  wave 
that  is  received  by  an  observer  moving  relative  to  the  source  of  the  waves.  The  relation 
between  the  actual  frequency  (/0)  and  observed  frequency  (/)  for  a  moving  target  is 
given  as 


/  =  /o 


f  v  \ 

1  +  -2- 

V  V  J 


(3.7) 


where  v  is  the  speed  of  the  waves  in  the  medium  and  va  is  the  speed  of  the  observer 
(positive  towards  the  source  and  negative  if  moving  away).  This  yields  an  increase  in  the 
observed  frequency  if  the  observer  is  moving  towards  the  source.  If  the  propagation 
medium  is  a  LHM,  this  effect  is  reversed  due  to  the  anti -parallel  phase  velocity  [6]. 
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Cherenkov  radiation  is  EM  radiation  emitted  when  a  charged  particle  moves  at  a 
speed  greater  than  that  of  light  in  the  medium.  The  radiation  propagates  in  the  forward 
direction  in  a  RHM.  In  a  LHM,  radiation  flows  backward,  opposite  to  the  particle 
velocity  [6]. 

C.  APPLICATIONS  OF  DNG  MATERIALS 

Negative-index  materials  have  attracted  great  attention  recently.  It  may  be 
possible  to  build  new  devices  by  using  metamaterials,  like  new  antenna  configurations  or 
perfect  planar  lenses.  Some  possible  applications  of  metamaterials  are  [4]: 

1.  Highly  sensitive  antennas. 

2.  Phase  and  dispersion  compensation  phase  shifters. 

3.  Antenna  and  microwave-device  miniaturization. 

4.  Broadband  and  multi-band  antenna  design. 

5.  Planar  lenses  that  exceed  the  diffraction  limit. 

The  list  above  includes  only  a  few  examples  among  the  possible  applications.  As  the  new 
properties  of  metamaterials  are  discovered,  the  list  is  likely  to  be  extended. 
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IV.  ANALYSIS  OF  SCATERING  FROM  MULTI-LAYERED 

METAMATERIALS 


A.  INTRODUCTION 

In  this  chapter,  the  wave  matrix  approach  is  used  to  calculate  the  reflection  and 
transmission  coefficients  of  multilayered  structures,  including  both  metamaterial  and 
dielectric  layers.  These  equations  were  programmed  in  MATLAB  6.5  so  that  the  effects 
of  frequency,  angle  of  incidence  and  polarization  could  be  examined. 

B.  MATLAB  SOFTWARE  AND  SIMULATIONS 
1.  Lossless  Dielectric  Panel 

The  first  configuration  considered  is  a  single  dielectric  slab.  Figure  11  shows  a 
dielectric  panel  that  has  a  thickness  of  t  =  0. 1  m  and  a  relative  dielectric  constant  of 
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Figure  1 1 .  Dielectric  panel  with  thickness  t  and  permittivity  s  =  4s0 . 


Assuming  nonnal  incidence,  from  Equations  (2.13)  and  (2.14),  the  reflection  and 
transmission  coefficients  become  [2], 

r  _?7i-7o 

7i+7o 

2)  =  i + rt  =  - 
3 
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(4.2) 


and 


r  . 1  V4 

7o+7i  1  +  _L 
A 


3 


(4.3) 


(4.4) 


It  is  also  possible  to  write  the  wave  matrix  for  the  panel  by  using  Equation  (2.19).  First, 
Figure  11  is  redrawn  to  show  the  incident  and  reflected  waves,  as  shown  in  Figure  12. 
Now,  the  wave  matrix  becomes 


^  A  if  e7'*"  r>- 


=  TT— 

A  I  r  ejm" 

n= 1  Ln  V1 


where  ®(;  =  yntn .  Since  the  thickness  of  the  last  layer  ( t2 )  extends  to  oo ,  it  is  possible  to 
arbitrarily  state  that  r3  =  0  and  ®2  =  0 .  Then,  is  obtained, 


>l_r  ejm 1  rA0iV  i  r2Y^3’ 

r  _  r  ^701  P-7®1  r  1  r 

Lm  J  v1  N  e  J\  2  1  y|_'3  _ 


(4.6) 


The  panel  transmission  and  reflection  coefficients  as  a  function  of  frequency  are 
plotted  for  the  range  of  500  MHz  to  2  GHz  in  Figure  13.  As  evident  in  Figure  13,  the  slab 
has  minimum  reflection  (maximum  transmission)  coefficients  at  0.75  and  1.5  GHz.  These 
are  the  frequencies  where  the  thickness  is  a  multiple  of  T/2 ,  where  X  is  wavelength  in 
the  dielectric. 

Figure  14  shows  the  transmission  and  reflection  coefficients  for  parallel 
polarization  at  various  angles  of  incidence  when  the  frequency  is  2  GHz.  The  slab  has  a 
minimum  reflection  coefficient  at  9j  =  64° ,  which  is  the  Breuster’s  angle. 

When  the  wave  is  perpendicularly  polarized,  the  slab  has  its  minimum  reflection 
coefficient  at  6j  =  0  .  In  this  case,  high-reflection  is  observed  at  grazing  angles  (Figure 
15). 
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Figure  12.  Incidence  and  traveling  waves  in  a  lossless  dielectric  panel  (After  Ref. 

[2])- 


Figure  13. 
a  =4 


Transmission  and  reflection  coefficients  of  a  lossless  dielectric  slab  with 
as  a  function  of  frequency  (nonnal  incidence,  parallel  polarization). 
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Figure  14.  Transmission  and  reflection  coefficients  of  a  lossless  dielectric  slab  as  a 
function  of  incidence  angle  (  /  =  2  GHz,  parallel  polarization). 
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Figure  15.  Transmission  and  reflection  coefficients  of  a  lossless  dielectric  slab  as  a 
function  of  incidence  angle  (/  =  2  GHz,  perpendicular  polarization). 
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2.  Radar  Absorbing  Material 

There  are  several  methods  of  reducing  the  RCS  (radar  cross  section)  of  targets. 
Radar  absorbing  materials  (RAM)  is  one  of  them.  The  next  case  examined  is  a 
combination  of  dielectric  and  magnetic  material  used  in  RAM  applications.  The  proper 
thickness,  material  configuration  and  electromagnetic  parameter  values  of  a  single  layer 
RAM  over  a  perfect  conductor  can  produce  zero  specular  reflection.  There  are  two 
different  approaches  to  achieve  zero  specular  reflection  [9].  The  first  is  the  matched 
characteristic  impedance  method,  in  which  the  characteristic  impedance  of  the  material  is 
made  equal  to  the  characteristic  impedance  of  the  free  space.  The  second  is  the  matched 
wave  impedance  method,  in  which  the  wave  impedance  of  the  material  is  made  equal  to 
that  of  free  space. 

Figure  16  shows  a  set  of  universal  curves  for  RAM  design  based  on  the  matched 
wave  impedance  method.  There  are  six  independent  parameters  incorporated  in  the 
universal  curves:  t.  A,  //' ,  //",  c'r  and  £ " .  They  can  be  reduced  to  four  by  normalizing 
pennittivity  and  permeability  by  t/  A .  Define  the  following: 

x  =  (t//L)Re[/y]  =  (t/T)//,' 

y  =  0  M)  Iml a]  =  (*/  (4  7) 

a  =  (t/A)Re[£]  =  (t/A)£'r  1  ' 

h  =  (t/ A)  Im[A]  =  ( t/  A)s 

In  Figure  16,  red  lines  represent  values  of  a  and  blue  lines  represent  values  of  tan  Se, 
where  tan  8e  =  a/b  .  To  enter  the  chart,  a  and  the  loss  tangent  (tan  <5))  are  needed.  Then,  it 
is  possible  to  find  the  sets  of  x  and  y  for  zero  reflection  by  using  the  chart. 

As  an  example,  a  design  of  a  RAM  layer  at  10  GHz  with  a  maximum  thickness  of 
3  mm  is  proposed.  If  the  layer  is  comprised  of  a  dielectric  with  s'r  =  10  and  s"r  =  100  ,  it  is 

possible  to  use  the  universal  curves  to  determine  p!r  and  //," .  To  enter  the  chart,  it  is 
necessary  to  find  a  and  the  electric  loss  tangent.  By  using  Equation  (4.7), 
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(t/x)  Im  [ |j.  ] 


=  ^/,  =  0.03  m ->-^  =  0.1 

(4.8) 

a  =  -U;  =  (o.ixio)  =  1 

A 

(4.9) 

_  ^  /  s" 

tan£  =a/u  =  —  =  10. 
e  /b  r’ 

(4.10) 

r 


By  going  to  the  intersection  of  curves  for  these  values  of  a  and  tan  Se ,  it  is  then  possible 
to  determine  x  and  y: 


— //'  1  '  u  «10 

(4.11) 

(4.12) 
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Jur  =  Ju'-Ju"  =  io-jm. 


(4.13) 


This  is  a  matched  characteristic  impedance  case  where  sr  =  //, .  In  general,  the  upper  part 
of  the  chart  represents  the  matched  characteristic  impedance  absorbers. 

Now,  to  demonstrate  the  zero  specular  reflection,  the  reflection  and  transmission 
coefficients  of  the  RAM  designed  by  using  the  universal  curve  using  the  wave  matrix 
formulation  will  be  plotted.  Figure  17  shows  reflection  and  transmission  coefficients  of 
the  RAM  as  a  function  of  angle  of  incidence  for  a  parallel  polarized  wave.  From  normal 
incidence  up  to  6j  =  35° ,  the  reflection  coefficient  is  smaller  than  0. 1 . 


Figure  17.  Transmission  and  reflection  coefficients  of  a  RAM  layer  as  a  function  of 

angle  of  incidence  (parallel  polarization). 


Figure  18  shows  the  reflection  and  transmission  coefficients  of  the  RAM  layer  for 
a  perpendicularly  polarized  wave.  As  evident  in  the  figure,  coefficients  are  polarization 
independent  because  /ur  =  sr  in  the  material. 
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1 


Figure  18.  Transmission  and  reflection  coefficients  of  a  RAM  layer  as  a  function  of 
angle  of  incidence  (perpendicular  polarization). 

3.  Propagation  in  Metamaterial  Slabs 

The  wave  matrix  formulation  for  multilayered  structures  is  given  by  Equations 
(2.19)  through  (2.26).  These  equations  are  based  on  right-handed  (RH)  materials.  The 
expressions  for  reflection  and  transmission  coefficients  are  identical  for  dielectric  and 
metamaterials.  On  the  other  hand,  the  phase  ®  in  Equation  (2.26)  is  preceded  by  a 
positive  sign  for  a  metamaterial  slab  and  by  a  negative  sign  for  a  dielectric  slab  [10]. 

First,  the  propagation  of  a  parallel-polarized  wave  in  a  single  metamaterial  slab 
will  be  studied.  To  verify  that  transmission  and  reflection  coefficients  computed  by  the 
software  are  correct,  the  slab  in  Figure  1 1  will  be  used,  but  this  time,  with  a  dielectric 
constant  4.  Figure  19  shows  transmission  and  reflection  coefficients  of  the 

metamaterial  slab  as  a  function  of  frequency.  The  frequency  range  covers  500  MHz  to 
750  MHz.  By  comparing  it  with  Figure  13  (dielectric  panel),  over  the  same  frequency 
range,  note  that  the  coefficients  are  the  same. 
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Figure  19.  Transmission  and  reflection  coefficients  of  a  metamaterial  slab  as  a 
function  of  frequency  (normal  incidence). 

a.  Antireflection  Coatings 

Next,  a  structure  consisting  of  LFI  and  RFI  slabs  in  air  (iV=4)  [10]  will  be 
studied.  Figure  20  shows  a  pair  of  metamaterial  and  dielectric  slabs  having  the  same 
width  ( t  =  0. 1  m),  and  opposite  permittivities  ( s  =  4,  er  =  -4  )  embedded  in  air.  In  this 

case,  rr  =0  and  tt  =  1  for  any  frequency  and  any  angle  of  incidence. 


In  Figure  21,  the  transmission  and  reflection  coefficients  of  the  structure 
are  plotted  as  a  function  of  frequency  for  a  parallel  polarized  wave,  and  Figure  22  shows 
coefficients  as  a  function  of  angle  of  incidence.  At  1  GHz,  the  advantage  of  structures  of 
alternating  slabs  is  that  the  total  reflection  and  transmission  coefficients  are  independent 
of  frequency  and  angle  of  incidence.  Adding  new  dielectric  and  metamaterial  pairs, 
which  are  identical  to  the  first  pair,  will  not  change  the  result.  Figure  23  shows  the 
reflection  and  transmission  coefficients  for  a  perpendicularly  polarized  wave.  As  evident 
in  the  figure,  the  structure  is  also  polarization  independent. 
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Figure  20.  Dielectric  and  metamaterial  slabs  embedded  in  air  (antireflection 

coatings). 
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Figure  2 1 .  Reflection  and  transmission  coefficients  of  a  pair  of  dielectric  and 
metamaterial  slabs  embedded  in  air  as  a  function  of  frequency  (normal  incidence,  parallel 

polarization). 


28 


1.4 


1.2 


|5  0.8 

% 

0  0.6 


0.4 


0.2 


0 


Reflection  Coefficient 
Transmission  Coefficient 


10  20  30  40  50 

Angle  (Degree) 


60 


70 


80 


90 


Figure  22.  Reflection  and  transmission  coefficients  of  a  pair  of  dielectric  and 
metamaterial  slabs  embedded  in  air  as  a  function  of  angle  of  incidence  ( /  =  1  GHz). 
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Figure  23.  Reflection  and  transmission  coefficients  of  a  pair  of  dielectric  and 
metamaterial  slabs  embedded  in  air  as  a  function  of  frequency  (normal  incidence, 

perpendicular  polarization). 
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The  phases  of  the  second  and  third  layers  are  the  same  but  have  opposite 
signs.  From  Equation  (2.26) 

$2=y^,A-sin2^  (4-14) 

and 

$3  “sin23  (4-15) 

where  t2  =  t3,  jur  =  /ur  and  sr^  =  -sr  .  Therefore,  theoretically  the  total  phase  of  the 

structure  must  be  zero.  Figure  24  shows  the  phase  of  the  total  transmission  coefficient  of 
the  structure  as  a  function  of  frequency.  As  evident  in  the  figure,  the  phase  of  the 
transmission  coefficient  is  zero  except  for  numerical  errors. 


Figure  24.  Phase  of  the  total  transmission  coefficient  for  a  normally  incident  parallel 

polarized  wave. 
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b.  High-Reflection  Coatings 

In  order  to  achieve  high-reflection,  metamaterial  and  dielectric  slabs 
should  be  chosen,  as  given  in  reference  [10] 

£,.  <£„,£„  >£,..£,.  <£,  (4.16) 

'A  r3  '3  r2  r2  r\  x  y 

and  02  =03  =n/2  at  the  central  frequency.  Figure  25  shows  a  pair  of  dielectric  and 

metamaterial  slabs  with  different  thicknesses  that  are  embedded  in  air.  The  permittivities 
are  chosen  to  maximize  the  reflection  coefficient.  In  this  case,  for  a  wave  that  is 
normally  incidence,  the  value  of  reflection  coefficient  varies  between  0.7  and  0.88 
depending  on  the  frequency  (Figure  26).  At  /  =  1.25  GHz,  very  high  reflection  values 
are  achieved  for  most  incidence  angles  (Figure  27). 


Figure  25. 
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Dielectric  and  metamaterial  slabs  with  different  thicknesses  embedded  in 
air  (high-reflection  coatings). 
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Figure  26.  Reflection  and  transmission  coefficients  of  a  pair  of  dielectric  and 
metamaterial  slabs  with  different  thicknesses  embedded  in  air  as  a  function  of  frequency 

(normal  incidence). 
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Figure  27.  Reflection  and  transmission  coefficients  of  a  pair  of  dielectric  and 
metamaterial  slabs  with  different  thicknesses  embedded  in  air  as  a  function  of  angle  of 

incidence  (/  =  2  GHz). 
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4.  Pairing  an  Epsilon-Negative  Slab  with  a  Mu-Negative  Slab  for  Zero 
Reflection 

An  Epsilon-negative  (ENG)  medium  has  negative  real  permittivity  and  positive 
real  permeability,  while  the  term  mu-negative  (MNG)  medium  is  used  for  the  medium, 
which  has  negative  real  permeability  and  positive  real  permittivity  [11].  It  is  possible  to 
achieve  zero  reflection  for  a  pair  of  ENG-MNG  slabs  with  specific  parameters  if  the 
incidence  wave  has  a  specific  angle  (Figure  28).  Zero  reflection  is  not  achievable  at  other 
angles  of  incidence  for  this  pair.  The  potential  application  of  a  matched  ENG-MNG  pair 
is  ideal  image  displacement  [11].  By  placing  an  object  in  front  of  this  pair,  an  observer  on 
the  other  side  of  the  pair  will  see  the  object  closer  by  the  amount  of  the  total  thickness  of 
the  pair. 


Figure  28.  A  pair  of  ENG-MNG  slabs  for  zero  reflection. 

The  parameters  for  the  ENG-MNG  slabs  are  given  in  Table  1.  The  propagation 
constant  is  given  in  Equation  (2.18).  Since  with  the  presence  of  loss  it  is  not  possible  to 
achieve  zero  reflection  due  to  the  mismatch  between  the  slabs,  it  is  assumed  that  ENG 
and  MNG  slabs  are  lossless.  For  this  matched  pair,  zero  reflection  is  achieved  only  for  an 
incidence  angle  of  45°  at /  =  3  GHz  (Figure  29). 
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Pennittivity  (<?) 

Permeability  (//) 

Thickness  ( t ) 

ENG 

-3^o 

2Ao 

2  n 

^  |  T ENG  1 

MNG 

60  s0 

-43.33 ju0 

2  n 

4-8 \/mng\ 

Table  1.  Parameters  for  a  matched  pair  of  ENG-MNG  slabs  (thin  slabs). 


The  reflectivity  depends  more  on  the  total  thickness  of  the  structure.  The  structure 
is  more  sensitive  to  the  variation  of  the  incidence  angle  for  a  thicker  pair  (Table  2). 
Figure  30  shows  the  reflection  coefficient  of  a  thicker  pair  as  a  function  of  the  angle  of 
the  incidence. 
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Figure  29.  Sensitivity  of  the  reflection  coefficient  to  the  variation  of  the  angle  of  the 

incidence  (thin  ENG-MNG  slabs). 
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By  comparing  Figure  29  and  Figure  30,  notice  that  the  reflection  coefficient 
remains  low  for  wider  range  of  incidence  angles  for  a  thin  pair.  The  reflection  coefficient 
increases  rapidly  as  the  angle  of  incidence  deviates  from  the  design  angle  ( 6i  =  45° )  in 
Figure  30. 


Permittivity  (<?) 

Permeability  ( //) 

Thickness  ( t ) 

ENG 

-3^o 

2Ao 

4  n 

5  / eng\ 

MNG 

60  s0 

-43.33 ju0 

2  n 

4-8  7mvg 

Table  2.  Parameters  for  a  matched  pair  of  ENG-MNG  slabs  (thick  slabs). 
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Figure  30.  Sensitivity  of  the  reflection  coefficient  to  the  variation  of  the  angle  of  the 

incidence  (thick  ENG-MNG  slabs). 
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The  potential  application  of  a  matched  ENG-MNG  slab  is  illustrated  in  Figure  3 1 . 
The  field  distribution  at  the  object  plane  can  be  expressed  by  Fourier  components  [11]. 
By  selecting  the  ENG  and  MNG  slab  parameters  ( sENG ,  /uh:NG ,  emng ,  juMNG  )  and  thicknesses 

( tENG ,  tMNG ),  the  Fourier  components  tunnel  through  this  matched  pair,  and  the  pair 

becomes  transparent  [11].  This  means  that  the  propagating  waves  will  go  through  the  slab 
and  show  up  at  the  exit  side  with  the  same  magnitude  and  phase.  An  observer  on  the  exit 
side  of  the  pair  will  see  a  virtual  image  of  the  point  source  as  though  it  was  placed  closer 
by  the  amount  of  the  total  thickness  of  the  slabs  ( tENC  +tMNG  ).  This  effect  can  be  used  in 

image  reconstruction  and  resolution  enhancement.  The  planar  lenses  with  negative 
refraction  are  described  in  Chapter  III.  These  lenses  create  a  real  image  of  the  source, 
whereas  ENG-MNG  pair  forms  a  virtual  image. 
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◄ - - ► 

t ENG  t MNG 


Figure  3 1 .  Pair  of  ENG-MNG  layers  under  certain  conditions  may  provide  image 

displacement  and  virtual  image  reconstruction  (From  Ref.  [12]) 

5.  Wet  Radome  Effect 

Water  flow  due  to  the  rain  on  the  surface  of  a  flat  radome  may  cause  a  water  film 
up  to  1  mil  (0.025  mm)  thickness.  This  thin  film  cannot  be  seen  except  from  reflections  at 
wide  angles.  The  water  layer  thickness  ( t )  on  a  flat  panel  is  given  as  [13] 
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r 


(4.17) 


t  = 


V 


3  RLju 
gp  cos  6* 


v/3 

J 


where  L  is  the  plate  length,  R  is  the  rain  rate,  //  is  the  specific  viscosity  of  water,  g  is  the 
gravitational  acceleration,  p  is  the  density  of  water  and  6  is  the  inclination  angle  from 
vertical.  The  combination  of  water  layer  and  radome  forms  a  multi-layered  structure  that 

can  be  analyzed  using  wave  matrices.  Figure  32  shows  transmission  loss 

.  .2 

(-lOlog  r  dB)  as  a  function  of  water  layer  thickness  for  a  parallel  polarized  wave 

normally  incident  upon  a  planar  radome  and  water  layer  at  20  GHz.  The  radome  consists 
of  a  0.03  in.  thick  laminate  of  fiberglass  and  a  water  repellent  skin  of  Tedlar  [14].  At  20 
GHz,  the  refractive  index  of  the  radome  is  n  =  1.7-  y'0.008  and  the  refractive  index  of 
water  [15]  is  n  =  6.85-  /2.63  .  When  the  radome  is  present,  transmission  increases  due  to 
the  improved  impedance  match  between  radome  surface  and  water  layer. 
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Figure  32.  Transmission  through  radome  and  water  layer  as  a  function  of  water  layer 

thickness. 

In  Figure  33,  transmission  loss  as  a  function  of  angle  of  incidence  for  a  0.032  in. 
fiberglass  radome  (sr  =  4)  is  plotted  at  15.5  GHz.  Figure  34  shows  the  transmission  loss 
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Transmission  Loss  (dB)  Transmission  Loss  (dB) 


for  the  same  radome  with  a  0.002  in.  water  layer  (sw  =  49-y'34).  High  loss  is  observed 
in  the  case  of  a  water  layer. 
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Figure  33.  Transmission  loss  through  a  fiberglass  radome. 
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Figure  34. 


Transmission  loss  through  wet  radome. 
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C.  SUMMARY 

This  chapter  is  dedicated  to  simulating  the  propagation  in  structures  consisting  of 
dielectric  slabs,  metamaterials  or  both  metamaterial  and  dielectric  slabs.  Different 
methods  are  described  to  achieve  zero  reflection  or  high  reflection.  These  methods 
include  designing  radar  absorbing  materials,  pairing  dielectric  and  metamaterial  slabs  and 
pairing  an  epsilon-negative  slab  with  a  mu-negative  slab.  For  all  structures,  reflection  and 
transmission  coefficients  are  plotted  as  a  function  of  frequency  and  angle  of  incidence.  It 
has  been  found  that  all  these  structures  have  interesting  features.  Potential  applications  of 
these  structures  are  described,  such  as  reducing  RCS  of  a  target  or  image  reconstruction. 
Finally,  transmission  through  a  radome  in  rain  is  discussed. 
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V.  CONCLUSIONS 


A.  SUMMARY  AND  CONCLUSIONS 

Equations  for  transmission  and  reflection  coefficients  for  a  single  boundary 
between  two  dielectric  media  are  presented  and  the  wave  matrix  approach  for  layered 
media  described.  Unusual  properties  of  metamaterials  are  discussed  and  potential 
applications  of  DNG  materials  are  presented.  Propagation  characteristics  of  planar 
multilayered  structures  (including  both  metamaterial  and  conventional  dielectric  layers) 
are  examined  by  using  wave  matrix  approach. 

A  combination  of  dielectric  and  magnetic  material  used  in  RAM  applications  is 
examined.  The  matched  characteristic  impedance  method  is  used  to  achieve  zero  specular 
reflection,  and  it  has  been  found  that  transmission  and  reflection  coefficients  are 
independent  of  polarization.  The  conditions  for  zero  specular  reflection  have  been  used  to 
validate  the  MATLAB  code. 

Anti-reflection  coating  and  high-reflection  coating  applications  have  been 
proposed.  It  has  been  found  that,  for  a  structure  consisting  of  a  pair  of  metamaterial  and 
dielectric  slabs  embedded  in  air  having  the  same  width  and  opposite  permittivities,  the 
reflection  coefficient  is  zero  and  independent  of  frequency  and  angle  of  incidence.  These 
structures  could  be  used  as  antenna  radomes  [10].  It  has  been  also  found  that,  by 
choosing  a  pair  of  dielectric  and  metamaterial  slabs  with  highly  concentrated  refractive 
indices,  the  reflection  coefficient  is  maximized.  For  one  such  structure  modeled,  the 
reflection  coefficient  is  larger  than  0.7  for  a  large  frequency  range. 

A  pair  of  ENG-MNG  slabs  has  been  examined  for  zero  reflection.  The  variation 
of  the  reflection  coefficient  with  angle  of  incidence  is  less  sensitive  for  thinner  slabs  than 
for  the  thicker  ones.  For  a  thin  slab,  the  reflection  coefficient  remains  low  over  a  wider 
set  of  incidence  angles.  By  selecting  the  ENG  and  MNG  slab  parameters  and  thicknesses 
properly,  the  pair  becomes  transparent  to  the  wave  [11].  This  effect  can  be  used  in  image 
reconstruction  and  resolution  enhancement. 

Finally,  the  effects  of  rain  on  transmission  through  radomes  is  discussed.  It  has 
been  found  that  the  existence  of  rain  increases  the  transmission  loss  substantially. 
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B.  FUTURE  WORK 

In  this  thesis,  infinite  planar  structures  consisting  of  dielectric  and  metamaterial 
were  examined.  No  edge  effects  and  no  surface  or  traveling  waves  were  considered. 
Future  research  could  be  conducted  to  examine  the  behavior  of  surface  waves  and  edge 
diffraction.  It  would  also  be  interesting  to  explain  how  metamaterials  could  be  applied  to 
radar  cross  section  reduction  and  anechoic  chamber  design.  In  principle,  absorbing  DNG 
layers  could  also  be  used  for  tenninating  the  computational  domain  for  the  finite  element 
(FEM)  method  and  the  finite  difference  time  domain  (FDTD)  method. 
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APPENDIX.  MATLABCODE 


This  appendix  provides  the  MATLAB  code  used  to  calculate  the  reflection  and 
transmission  coefficients  of  multi-layered  structures  in  Chapter  IV.  The  code  is  written  in 
MATLAB  6.5  and  based  on  wave  matrices  model.  This  code  calculates  and  plots 
reflection  and  transmission  coefficients  as  a  function  of  frequency.  It  can  be  modified 
easily  to  calculate  the  coefficients  as  a  function  of  angle  of  incidence. 


%################################################ 

%SCATTERING  FROM  MULTI-LAYERED  METAMATERIALS 
%############################################### 

%This  program  calculates  and  plots  refraction  and  transmission  coefficients 
%of  multi-layered  structures  as  a  function  of  frequency  for  a  specific 
%angle  of  incidence 

%INPUTS 

%Defining  structure  properties 

clear; 

clc; 

N  =  input('Enter  the  number  of  layers=  '); 

fmin=input('Enter  Minimum  Frequency  (GHz)=  '); 

fmax=input('Enter  Maximum  Frequency  (GHz)=  '); 

th=input('Enter  Incidence  Angle  (Degree)=  '); 

p=input  ('Enter  Polarization  Type  Parallel  (1),  Perpendicular  (2)='); 

%Constants 

E0=8.85e-12;  %Pennittivity  of  free  space 
Mu0=12.5664e-7;  %Permeability  of  free  space 
c=l/sqrt(EO*MuO);  %Speed  of  light  (m/s) 

%Defining  layer  parameters 
for  layer=l:N  ; 
display(layer); 

Mu(layer)=input('Enter  Permeability(p)  = '); 

Ei(layer)=input('Enter  Permittivity(€)  = '); 
t(layer)=input('Enter  Layer  Length  (m)  =  '); 
sigma(layer)=input('Enter  Loss  (sigma)  =  '); 
end 

f=linspace(fmin,fmax);  %Frequency  matrix 
thr=th*pi/180;  %Angle  of  incidence  (Radian) 
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st=sin(thr); 

ct=cos(thr); 


for  i=l:N 
if  Ei(i)<0 
E(i)=abs(Ei(i)); 
else 

E(i)=Ei(i); 

end 

end 

%This  part  calculates  the  total  reflection  and  transmission  coefficients  of 
%the  structure 
for  k  =l:length(f); 
lambda=c/(f(k)*  le9); 

Beta  1  =2  *  (pi/lambda); 
total=l; 

w(k)=2*pi*f(k)*  le9; 

%Calculating  the  characteristic  impedance  of  the  layers 
for  i=l:N 

Er(i)=E(i)-(j*sigma(i)/w(k)/EO); 
if  p==l 

Z(i)=sqrt(Er(i)  *Mu(i)-stA2)/Er(i)/ ct; 
elseif  p==2 

Z(i)=Mu(i)*ct/sqrt(Er(i)*Mu(i)-stA2); 

else 

display  ('Polarization  type  must  be  1  or  2’); 
end 
end 

%Individual  reflection  and  transmission  coefficients  of  the  layers 

ro(  1  )=(Z(  1 )- 1  )/(Z(  1 )+ 1 ); 

tao(l)=ro(l)+l; 


for  i=2:N 

r°(i)=(Z(i)-Z(i- 1 ))/ (Z(i)+Z(i- 1 )); 
tao(i)=ro(i)+l; 
end 

%Generating  wave  matrices 
for  i  =1:(N-1) 
if  Ei(i)<0 
sign=- 1 ; 
else 
sign=l; 
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end 

Beta=Betal  *sign; 

A=exp(j*Beta*t(i)*sqrt(Er(i)-stA2)); 

Ac=exp(-j*Beta*t(i)*sqrt(Er(i)-stA2)); 

total=total*((l/tao(i))*[A,ro(i)*Ac;ro(i)*A,Ac]); 

end 

%Total  reflection  and  transmission  coefficients  of  the  structure 
total=l/tao(N)*total*[l  ro(N);  ro(N)  1]; 
rot(k)=total(2,  l)/total(l ,  1); 
taot(k)=  1  /total(  1 , 1 ) ; 
end 

%This  part  plots  the  coefficients  as  a  function  of  frequency 

plot(f,abs(rot),'r-',f,abs(taot),’— '); 

grid; 

xlabel('Frequency  (GHz)'); 
ylabel('Coefficient  Magnitude'); 

legend('Reflection  Coefficient', 'Transmission  Coefficient'); 
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